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Abstract 

The  properties  of  metal -metal  interfaces  are  of  interest  in  many  metallurgi- 
cal applications.  These  include  grain  boundaries,  crack  growth,  friction  and  bi- 
metallic adhesion.  Tne  present  work  is  a study  of  the  electronic  properties  of 
a simple  bimetallic  junction.  The  methodfe  employed  to  investigate  the  interface 
are  the  Green’s  function  technique  and  the  phase  shifts  method.  We  calculate  the 
Green's  function  of  a simple  junction  within  the  tight -binding  approximation.  The 
conditions  for  the  occurance  of  bound  states  is  deduced  from  the  poles  of  the 
Green's  function.  It  is  shown  that  there  are  three  types  of  wave  functions  asso- 
ciated with  the  interface.  The  first  extends  throughout  the  entire  crystal,  the 
second  extends  on  one  side  of  the  interface  only,  whereas  the  third  kind  is  localized 
near  the  interface.  Using  the  phase  shifts  method  we  derive  an  expression  for  the 
change  in  density  of  states  due  to  the  creation  of  the  interface.  From  this 
expression  we  derive  the  corresponding  single-particle  contribution  to  the  interface 
energy  and  the  interface  specific  heat. 
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I.  Introduction 


The  properties  of  metal -metal  interfaces,  are  of  interest  in  many  metallurgi- 
cal applications.  These  include,  for  example,  grain  boundaries,  crack  growth, 
friction  and  bimetallic  adhesion.  The  general  metal -metal  interface  is  a nuch 
more  complicated  problem  than  the  corresponding  metal -vacuum  interface.  Recently 
there  has  been,  however,  some  success  in  applying  certain  methods  developed  for 
surface  phenomena,  to  the  investigation  of  bimetallic  interfaces. 

The  electronic  properites  of  1-dimensional  bimetallic  junctions  were  dis- 
cussed by  several  authors  using  various  types  of  models.  Aerts^  studied  the 
electronic  structure  of  a 1-dimensional  Kronig-Penney  model,  in  the  limit  of  a 
6-like  potential.  In  this  work  he  established  the  possibility  of  the  existence  of 
bound  interface  states.  A tight-binding  approach,  which  was  used  successfully 

for  several  surface  problems,  was  also  applied  to  the  1 -dimensional  interface. 

(2) 

Davison  and  Cheng  ' investigated  the  electronic  properties  of  such  a system 

using  the  molecular-orbital  method.  The  nodel  they  use  associates  a single  s-type 

(3) 

orbital  with  each  atom.  A similar  model  was  studied  by  Allan  and  Larmoo  , using 

an  approximate, Gaussian,  density  of  states,  having  the  correct  second  moment. 

(4) 

Green's  function  formalism  was  used  by  Foo  and  Wong  y to  study  the  interface 
states  of  a 1-dimensional  sp-hybrid  junction. 

As  far  as  we  know,  the  study  of  three-dimensional  bimetallic  inte  vas 

carried  out  by  the  density-functional  formalism  only.  This  formalism,  dt.  sloped 
by  Hobenberg,  Kohn,  and  Sham,  ^ ^ was  applied  recently  to  bimetallic 
interfaces  by  Bennett  and  Duke,  ^ Ferrante  and 
and  Schuttler , and  by  Mehrotra , Pant  and  Das.^) 
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Whereas  the  density-functional  formalism  is  applicable  mainly  to  simple 
metals,  the  tight-binding  approximation  is  more  suitable  for  the  description 
of  transition  metals.  Since  the  purpose  of  the  present  work  is  to  investigate 
the  electronic  structure  of  a bimetallic  junction,  formed  from  two  transition 
metals,  we  shall  use  the  tight -binding  approach.  The  model  we  consider  is  a 
highly  simplified  one.  The  two  metals,  on  each  side  of  the  junction,  are  des- 
cribed by  s-type  tight-binding  Hamiltonians . The  interface  we  consider  is 
formed  by  bringing  together  two  semi- infinite,  simple  cubic  crystals,  and  creating 
bonds  between  the  atoms  on  the  two  sides  of  the  interface.  We  assume  that  the 
two  semi- infinite  crystals  have  the  same,  two-dimensional,  translation  symmetry 
parallel  to  the  interface.  The  electronic  properties  of  this  model  are  investi- 
gated by  using  the  Green's  function  method,  described  in  detail  by  Kalkstein  and 
Seven/13) 

The  details  of  the  tight -binding  model  are  outlined  in  Section  II.  Within 
this  model  we  allow  for  a change  in  the  self-consistent  potential  of  the  electrons 
near  the  interface.  The  diagonal  matrix  elements  of  the  Green's  function  are 
calculated  in  Section  III  by  considering  the  formation  of  the  interface  as  a 
perturbation  on  the  two  sani- infinite  crystals.  This  is  accomplished  by  the  use 
of  Dyson's  equation.  Using  the  expression  for  the  Green's  function  we  discuss 
the  electronic  structure  of  the  interface  in  Section  IV.  We  show  that  there  are 
three  types  of  wave  functions  associated  with  the  interface.  The  first  one 
extends  throughout  the  entire  system,  the  second  type  extends  on  one  side  of  the 
junction  only,  and  the  third,  associated  with  bound  states,  is  localized  near  the 
interface.  The  behavior  of  the  bound  states  as  a function  of  the  coupling  constant 
between  the  two  metals  is  also  discussed.  The  "phase  shifts"  method  is  applied 
in  Section  V to  determine  the  change  in  the  total  density  of  states  due  to  the 
creation  of  the  interface.  Applying  this  result  we  obtain  an  expression  for  the 
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single-particle  energy  contribution  to  the  interface  energy  and  the  corresponding 
contribution  to  the  electronic  specific  heat.  An  application  of  the  model  to  the 
interface  formed  between  two  transition  metals  of  the  same  series  is  presented 
in  Section  VI,  where  we  calculate  numerically  the  local  densities  of  states , 
the  change  in  density  of  states  and  the  interface  energy. 

II.  The  tfodel 

Consider  the  formation  of  a metal-metal  interface  by  bringing  together 
two,  semi- infinite,  metallic  crystals.  As  soon  as  a contact  is  formed,  electrons 
will  flow  from  the  metal  having  the  higher  Fermi  energy  to  the  one  having  the 
lower  energy.  This  flow  of  electrons  stops  when  the  potential  energy  difference 
between  the  two  sides  of  the  junction,  which  is  created  by  the  dipole  layer 
produced  at  the  interface,  is  equal  in  magnitude  and  opposite  in  sign  to  the  difference 
between  the  two  Fermi  levels.  Thus,  the  Fermi  levels  of  the  two  metals  are  aligned, 
due  to  the  interface  dipole-layer,  when  the  bimetallic  junction  is  formed.  If 
-Av  and  +Av  are  the  electrostatic  potential  created  by  the  dipole  layer  on 
the  right  and  the  left  hand  sides  far  away  from  the  interface,  then  we  have: 


eAv  = % (Efa  - Efb)  (2.1) 

ci  b 

where  -e  is  the  charge  of  the  electron  and  E^  and  Ej  are  the  original 
Fermi  levels  of  metals  a and  b , located  to  the  right  and  to  the  left  of 
the  junction,  respectively.  The  cannon  Fermi  level,  after  the  junction  is  formed, 
is  given  by: 


Ef  = Efa  - eAv 
= Efb  + eAv 

= h (Efa  + Efb)  (2.2) 


In  order  to  describe  the  electronic  properties  of  the  tvro  bulk  metals,  the 
two  metal -vacuum  interfaces,  and  the  metal -metal  interface,  we  apply  the  tight- 
binding  approximation  for  the  various  Hamiltonians.  For  simplicity  we  associate 
one  s-type  orbital  with  each  lattice  site.  It  is  also  assumed  that  each  orbital 
has  a g-fold  degeneracy  in  order  to  partially  account  for  the  10-fold  degeneracy 
of  the  d orbitals  in  transition  metals.  Therefore,  the  Hamiltonians  of  the 
two  bulk  metals,  Ha  and  H^  , are  given  respectively  by 


H 

a 


_■  u 

h v.'ij 

1 , J 


t 

b .b  . 
i J 


(2.3) 


where  a., a.  and  b.,  b.  are  the  creation  and  destruction  operators  of  the 
i j i j 

Warmier  type  orbitals,  associated  with  metals  a and  b(  localized  near  sites  i 

and  j . The  prime  on  the  summation  signs  denotes  a summation  over  nearest  neighbors 

o b 

only.  The  matrix  elements  ty  and  ty  are  given  by 


a 

t . . = 


-ij 


ci!Ha,j>  = 


- 


l«b'J>  * 


E if  i = j 
oa  J 

E^a  or  E*a  if  i and  j are  nearest  neighbors 

0 otherwise 

5 . if  i = j 
ob  * 

E^  if  i and  j are  nearest  neighbors 

0 otherwise 


(2.4a) 


(2.4b) 
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To  cleave  the  crystal  along  a given  low-order  crys talographic  plane,  we  have  tc 
break  the  bonds  between  two  adjacent  planes,  parallel  to  the  corresponding  direc- 
tion. We  neglect  any  geometrical  reconstruction  of  the  crystal  due  to  the  forma- 

sl  b 

tion  of  the  surface  and  assume  that  the  transfer  integrals  ty  and  t^  with 
i f j , have  the  same  value  as  in  the  bulk  systems,  provided  that  both  sites 
i and  j are  occupied.  In  addition,  we  assume  that  the  effect  of  the  surface 
on  the  redistribution  of  the  electronic  charge  near  the  surface  can  be  described  as 
a change  in  the  self-consistent  potential  of  the  electrons  near  the  first  surface 
layer  only.  Therefore,  the  diagonal  matrix  elements  of  the  Hamiltonian  for  sites 
located  on  the  surface  plane  will  be  different  from  the  corresponding  bulk  value. 
Thus,  for  the  metal-vacuum  systems  we  have 


<i|H;|i>  = E 
*-i|H^|  i>  = E 


oa 


ob 


+ U 

a 


+ cb 


6. 


-1 


(2.5a) 

(2.5b) 


where  Ua  and  denote  the  change  in  the  self-consistent  potential  of  the 
electrons  near  the  surface.  In  the  above  expressions  we  assumed  that  i = 0 and 
i = -1  denote  surface  sites  of  metal  a and  b , respectively. 

In  order  to  form  the  metal -metal  interface  we  start  from  two  non-interacting 
semi- infinite  metal -vacuum  interfaces  (the  free  surfaces) , having  the  same  trans- 
lational synmetry  parallel  to  the  surface,  and  introduce  a coupling  between  the 
two  surface  layers.  The  Hamiltonian  of  the  non-interacting  system  is  given  by 


H°  = 


ct 

i 


c. 

J 


(2.6a) 


where 


t. . = 


<i|H'|j>  if  both  i and  j are  a sites. 

a 


ij  ~ j<i|H^)|j>  if  both  i and  j are  b sites. 


o otherwise 


c c+  = 
i’  i 


ai  a/  ^ i is  311  a site 


4. 

b^,  b/  if  i is  a b site . 


The  interaction  between  the  two  surface  layers  is  characterized  by  a transfer 
matrix  element 

•>'< 

< R JH|R  > = <R  |H|R  ,>  <2-7 

— 1 — o — o'  1 — l 

where  R and  R Q denote  two  nearest-neighbors  sites  located  on  the  two  sides 
of  the  interface.  We  also  note  that  the  dipole  layer  created  at  the  interface 
will  change  the  self-consistent  potentials  on  the  two  sides  of  the  junction. 

If  we  assume  that  charge  redistribution  is  confined  to  the  near  vicinity  of  the 
interface  we  can  write 


<-R.  |H  Ir  > = e - eAV  + U 6.  , if  i is  an  a site 

1 1 oa  a i,o 


<Sj ! H | Rj > = Eob  + eAV  + Ub  5^  ^ if  j is  a b site 


where  i = 0 and  j = -1  denote  interface  sites  of  metals  a and  b , respec- 
tively. In  the  above  expressions  and  denote  the  change  in  the  self- 
consistent  potential  of  electrons  near  the  interface  relative  to  the  bulk. 

Ill . The  Bimetallic  Interface  Green's  Function 

As  we  already  noted  before,  we  assume  that  the  two  crystals  forming  the  inter- 
face have  the  same  translational  symmetry  parallel  to  the  interface  plane.  As 
a result,  the  wave  vector  parallel  to  the  interface  k(|  is  a good  quantum  number. 

In  the  following  we  assume  that  the  interface  is  in  the  zy  plane,  and  that 
metals  a and  b are  to  right  and  to  the  left  of  the  interface  respectively. 

The  integers  m,  n,  ...  will  be  used  to  label  the  various  planes  parallel  to 
the  interface.  Metal  a is  assumed  to  occupy  the  planes  m _>  0, whereas  m £ -1 

(13) 

planes  are  occupied  by  b atoms.  Following  Kalkstein  and  Soven  ' , we  denote 
by  xn  the  translation  vector  parallel  to  the  interface  which  brings  the  atoms 
in  the  n-tb  - lane  to  coincide  with  the  transverse  positions  of  the  atoms  on  one  plane 
n = 0 neral  coordinate  of  an  atomic  site  on  the  n-th  plane  is  thus  given 

by 


R = R i + R ..  + 1 
~n  "n  n 


(3.1) 


where  Ri  is  the  distance  between  the  planes  o and  n , and  R,,  is  a general 

translation  vector  parallel  to  the  interface.  Using  the  localized  Wannier  orbitals 

|R  i +R,,  + xn>  we  define  the  mixed  Bloch-Wannier  representation  by  the  following 

(14) 

two-dimensional  Bloch  sunr 


k,,,n>  = (N'„) 


-*s 


I 

R.i 


R„ 


+ T " 

~n 


e%- 


CR,, 


In> 


(3.2) 
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where  N„  is  the  lumber  of  atoms  in  the  plane  parallel  to  the  interface.  These 
functions  are  localized  near  the  n-th  plane.  Since  the  translational  synmetry 
parallel  to  the  interface  is  preserved,  the  eigenfunctions  of  the  interface 
Hamiltonian  can  be  expanded  in  terms  of  the  Bloch-Wannier  functions  (3.2)  as 
follows 


|E  k„>  = £ C~\  |k„,  n> 

j n n£-. 


(3.3) 


where  j stands  for  any  quantum  number  needed  to  specify  the  eigenstates  in 


addition  to  kM 


Using  the  expression  (3.3)  we  can  define  various  densities  of  states  for 
the  system  with  the  interface.  The  local  density  of  states,  near  the  n-th  plane, 
with  a definite  value  of  k„  is  given  in  terms  of  the  coefficients  (P^  by 


0„  - r |CJ  1 t(E-E  > 


The  local  density  of  states  (per  atom)  near  the  n-th  plane  is 


(3.4) 


(E)  N ^ Pn  (E,k„)  N 
n , n 1 " i u 

K 1 1 J » fill 


C , 6 (E-E  , ) 


(3.5) 


The  son  over  k in  this  expression  goes  over  all  values  of  k M in  the  two- 
dimensional  Brillouin  zone  defined  by  the  crystals  geometry  parallel  to  the  inter- 
face. Summing  the  expression  (3.5)  over  the  various  crystal  planes  we  obtain 
the  total  density  of  states.  The  density  of  states  per  atom  is  thus  given  by 


1 6(E  - Ejk„> 


(3.6) 


p(E) 


(E) 


C\ 

nk„ 


where  N 1 is  the  nunber  of  atomic  planes  parallel  to  the  interface  and  N is  the 
total  number  of  atoms  given  by  N = N„  N i, 

In  order  to  derive  the  various  densities  of  states  of  the  bimetallic  inter- 
face we  shall  use  the  Green's  function  technique.  Let  H and  H°  be  the  Hamil- 
tonians of  the  interfaced  crystal  and  of  the  free  metal -vacuum  interfaces  respec- 
tively. The  metal-metal  interface  Green's  operator  G is  defined  by  the  following 
equation 


(E  - ifi  - H)  g = 1 


(3.7) 


where  E is  the  energy  and  6 is  a positive  infinitesimal.  The  surface  Green's 
operator  G°  is  related  to  the  Hamiltonian  H°  by  a similar  equation.  The 
Green's  operator  of  the  interface  is  related  to  the  surface  Green's  operator  via 
Dyson's  equation: 


G = G°  + G°VG 


(3.8) 


where  V is  the  perturbation  necessary  to  create  the  metal -metal  interface  from 
the  free  surfaces  of  metals  a and  b , i.  e. 


I 


We  note  that  because  of  the  translational  symmetry  parallel  to  the  interface, 

G,  G°  and  V will  all  be  diagonal  in  the  wave  vector  index  k,,  , in  the  Bloch- 
Wannier  representation  (3.2).  Omitting  the  corresponding  6-function, 

6(k,,  - k,.’),  we  use  the  notations  G (m,n;  k,.)  , G°  (m,n;  k,,)  and  V (m,n;  k,,) 
for  the  matrix  elements  of  G,  G°,  and  V in  the  Bloch-Warmier  representation. 

To  further  simplify  the  notation,  we  shall  generally  omit  the  explicit  k,,  depen- 
dence in  these  expressions. 

The  various  densities  of  states  defined  earlier  by  Eqs.  (3.4)  - (3.6),  are 
simply  related  to  the  imaginary  part  of  the  diagonal  matrix  elements  G (m,m;  k,,) 


,(E)  = ImTr  G = 4 In  E G (bi.i-.i;  k„) 


3 I'E')  = — in'  £ G(n,n;  k„) 
n v ttN„  v 


(E,  k„)  = p ImG  (n,n;  k„) 


(3.10) 


(3.11) 


(3.12) 


We  turn  now  to  the  evaluation  of  the  interface  Green's  function  from  Dyson's 
equation  (3.8).  This  operator  equation  is  reduced  to  the  following  algebraic 
equation  in  the  Bloch-Warmier  representation: 


G(m,n)  =G°(m,n)  + £ G°  (m,  4)V  (■£  ,r  )G  (r  ,n) 

i,r 


(3.13) 


From  the  discussion  of  the  model  given  in  Section  II,  it  is  obvious  that  the 
perturbation  potential  V has  off-diagonal  matrix  elements,  which  couple  the 
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interface  planes  n = 0 and  n = -1  . In  addition,  on  each  side  of  the  inter- 
face V has  two  types  of  diagonal  matrix  elements.  The  first  one,  due  to  the 
interface  dipole  layer,  is  given  by  +eAv  . This  is  just  the  perturbation 
necessary  to  align  the  Fermi  levels  on  the  two  sides  of  the  junction,  Eq.  (2.1). 

The  second  type  of  diagonal  matrix  elements  is  due  to  the  change  in  the  self- 
consistent  potential  of  the  electrons  near  the  interface.  This  is  given  by 
Ua  - Ua  and  where  U , Ua  and  U^,  are  the  self-consistent  potentials, 

relative  to  the  bulk,  of  interface  and  surface  electrons  of  metals  a and  b , 
respectively. 

As  we  shall  see  later  on,  it  is  possible  to  determine  the  surface  Green's 
functions  for  an  arbitrary  value  of  the  surface  self-consistent  potentials  Ua  and 
. It  is,  therefore,  simpler  to  start  from  a fictitious,  intermediate  surface 
problem  where  the  surface  self -consistent  potentials  Ua  and  have  already 
the  correct  interface  values,  Ua  and  , respectively.  Using  this  system 
as  our  starting  point,  the  perturbation  potential  necessary  to  create  the  metal- 
metal  interface  is  simpler  than  in  the  original  problem.  It  consists  of  diagonal 
terms  ±eAv  , which  align  the  Fermi  levels  of  the  two  metals  and  of  an  off-diagonal 
term,  which  couple  the  two  semi- infinite  crystals.  Thus,  the  only  non-vanishing 
matrix  elements  of  the  perturbation  are  given  by 


V(m,m)  - 


-eAv  for  m ' 0 
+eAv  for  m s -1 


(3.14) 


V(-1,0)  = V (0,-1)*  = 9(k.,)  e1<t>(-"^ 


(3.15) 


where  0 and  <p  are  derived  from  the  relation 

V(-1,0)  - L'  <R  f Vf R > ‘ + l) 

R 1 -O 

-o 

In  this  expression,  t is  a translation  vector  parallel  to  the  interface 
which  brings  the  transverse  atomic  sites  of  the  two  interface  planes  to  coincide. 

In  Eq.  (3.15)  the  phase  <J>  was  chosen  in  such  a way  that  0>O. 

We  start  now  from  the  intermediate  surface  problem  discussed  above,  whose 
Green's  function  we  denote  by  G°  , and  apply  the  perturbation  in  two  steps.  First 
we  align  the  Fermi  levels  of  the  two  metals  by  applying  a constant  electrostatic 
potentials  +Av  on  metal  a and  -AV  on  metal  b . The  only  effect  of  applying 
these  constant  potentials  on  the  surface  Green's  functions  is  a shift  in  the 
corresponding  energies.  Explicitly  we  have 

G°  (E)  = G°  (E-hjAv) 

G°  ( E ) = G°  (E-eAv) 

where  G denotes  the  value  of  the  Green's  functions  in  the  intermediate  problem, 
after  the  application  of  the  potentials  ±Av  . To  simplify  the  notation  we  drop 
the  wiggle  from  these  Green's  functions  remembering  that  tine  energies  have  to  be 
shifted  according  to  (3.17). 

At  this  stage  we  turn  on  the  second  part  of  the  perturbation  i.  e.  the  off- 
_aronal  coupling  between  the  two  metals.  Due  to  the  localized  nature  of  this 
perturbation,  the  Dyson's  equation  (3.13)  is  greatly  simplified,  and  we  obtain 

G(m,n)  =G°(th,n)  +G°(m,-1)  V(-1,0)  G(0,n)  +G°(m,0)  V (0,  -1)G (-1  ,n)  (3-18) 


(3.17a) 

(3.17b) 
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Using  the  fact  that  the  surface  Green's  function  G° (m,n)  vanishes  if  m and 
n refer  to  planes  on  opposite  sides  of  the  interface,  we  can  easily  solve  equa- 
tion (3.18)  for  the  perturbed  Green's  function.  In  this  way  we  get 


C(m,n)  =G°(m,n)  +G°(m,0)  [v(0,-l)G°(-l,n)  + | V(0, -I)  | 2G°  (-1, -1) 
G° (0,n)]  [l  - | V (0, -1) i 2 G°(0>0)G°(-ll-l)J 


-1 


(3.19a) 

for  m>0 
and 

G(m,n)  =G°(m,n)  + G°  ( m,-l)  [ V(-1,0)  G°(0,n)  + ! V (0, -1) | 2G° (0,0)  x 
G°(-l,n)J  [l-  | V (0, -1) i 2 G°(0,0)  G°  (-1 , — 1 ) J 

(3.19b) 

or  m _<  -1. 

In  particular,  the  diagonal  matrix  elements  of  the  interface  Green's  function  are 
given  by 

-1 


G(m,m)  «G°(m,m)  +B2G°(-1,-1)  G°  (m,  0)G°  (0  ,m)  | 1-bV  (0 ,0)G°  (-1, -1  ) 


for  m >_  0 
and 


(3.20a) 


G(m,m)  = G (m 


,„>  +82G»(0.o)  G°(»,.l)  G° | 1-62  G°(0.0)GO(-l.-l)  | 


-I 


for  m < -1. 


(3.20b) 
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In  the  above  expressions  g is  related  to  the  off-diagonal  coupling  between  the 
two  metals  through  Eq.  (3.15),  i.  e. 

B - | V(0,-1) | 

The  diagonal  matrix  elements  of  the  surface  Green's  function  were  calculated 
by  Kalkstein  and  Seven . Using  the  same  method  one  can  also  derive  the 
general  matrix  elements  of  this  Green's  function.  In  Appendix  A we  show  that 
these  matrix  elements  are  given  by 


*m-n  I ;v 


i^+(u>  -2U  ) 
d a a 

iu-(u)  -2U  ) 
a a a 


for  m,  n>  0 


0 -1  -t(m-n)6.  I /VLlJb 

G (tn.n)  = e b j \ 2Tb 


n I -2 


m • 


WV20b) 


n*<V2V 


for  m,n  < -1 


G (ro,n)  - 0 otherwise. 


(3. 


In  these  expressions  y and  oj  are  defined  by 

- 

(4T2-u>2)^  for  uj2  £.  4T2 
isign(uj)  (uj2-4T2)^  for  w2  > 4T2 

and 

OJ  = E - W(k„) 


(3.21c) 


where  W,  T,  and  9 are  related  to  the  matrix  elements  of  the  bulk  Hamiltonians 
of  the  two  metals, in  the  Bloch-Warmier  representation,  as  follows 


Wa(b)  (-n)  = n-"'tla(b)  Ink,, 

Ta  (b ) m)  ei^a (b)  = <nk„|Ha 


a (b)  ! n+1  ,klr> 


In  order  to  apply  the  expressions  (3.21)  to  the  solution  of  the  interface 
problem,  we  replace  first  UQ  and  by  the  corresponding  interface  values 
Ua  and  , respectively.  In  addition  we  have  to  shift  the  energies  according 
to  Eq.  (3.17)  in  order  to  align  the  Fermi  levels  of  the  two  metals.  This  can 
be  achieved  by  redefining  the  w-s  in  the  following  way 


xa  ~ E - wa(k.i)  + eAv 
^b  = E " Wb^-")  " eAv 


(3.22a) 

(3.22b) 
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Substituting  the  explicit  expression  (3.21)  for  the  surface  Green's  function 
into  Eq.  (3.20),  we  obtain  that  the  diagonal  matrix  elements  of  the  interface 
Green's  function  are  given  by 


G(m,m)  - — + 


. +L(w 
a a a 


l /_V^a\ 

-2U  ) \ 2T  / 


2m  / 


i^a+  (V  2Ua) 


2 

89  i 


(3.23a) 


4b2+  J_na+  | 


for  m > 0 , and 


2 m -2 


G(m,m)  = — + 


■+!'(«*. 


— (=£H  ! 

U)  -2UU)  b 


^+<V2Ub> 


892i 


J 


402  + [ua+  i(UJa-2Ua)]_[^+  t(»b-2Ub)  j 


(3.23b) 


for  m ^ -1. 

We  recall  tliat  wa  and  in  these  equations  are  defined  by  the  energy-shifted 

expressions  (3.22a)  - (3.22b),  respectively. 

It  is  straightforward  to  confirm  that  in  the  limit  where  the  two  metals 
a and  b are  the  same,  the  expression  (3.23)  for  the  interface  Green's  function 
reduces  to  the  corresponding  bulk  expression. 

IV.  The  Electronic  Structure  of  the  Interface 


In  this  Section  we  apply  the  Green's  function  derived  in  the  preceding  section 
in  order  to  investigate  the  electronic  structure  of  the  interface.  From  the  expli- 
cit expression  (3.23)  it  follows  that  the  diagonal  matrix  elements  of  the  inter- 
face Green's  function  have  a non- vanishing,  continuous  imaginary  part  for  energies 
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which  are  either  insid*,  the  band  of  metal  a or  inside  the  band  of  iretal  b 
(shifted,  of  course,  according  to  (3.17)).  Therefore,  as  follows  from  Eq.  (3.12) 
the  band  width  of  the  crystal  with  the  interface  is  the  union  of  the  shifted 
bands  of  the  two  separate  crystals.  However,  the  wave  functions  of  the  combined 
crystal  can  be  classified  into  three  distinct  classes  according  to  their  different 
localization  properties.  To  facilitate  our  further  discussion  we  define  the  k,, 
sub-bands  of  the  two  metals  as  the  band  structure  obtained  by  the  intersection  of 
Ea(k)  and  E^k)  with  V.  = constant,  where  Ea(k)  and  E,(k)  are  the 

single-particle  energy  s>  >e  two  bulk  crystals. 

The  first  type  of  states  nave  /rave  functions  which  extend  throughout  the 
entire  crystal.  As  can  be  seen  from  Eqs.  (3.4)  and  (3.23),  this  kind  of  behavior 
is  associated  with  states  whose  energy  lies  in  the  k,,  sub-band  of  the  two  metals, 
i.e.  their  energy  and  wave  vector  satisfy  the  relations 

Ie  - Wa  (k„)  + eAv  | <2T  (k„) 

d a 

|e  - Wb(k„)  - eAv|  <2Tb(k„)  (41) 

The  second  type  of  states  have  wave  functions  which  extend  to  infinity  on 
only  one  side  of  the  interface,  and  which  decay  exponentially  with  the  distance 
frcm  the  interface  on  the  other  side.  Using  the  expression  (3.23)  for  the  Green's 
function  it  is  not  hard  to  see  that  this  behavior  occurs  when  the  energy  lies  in 
the  k,.  sub-band  of  one  metal  but  outside  the  corresponding  sub-band  of  the  other. 
Thus,  for  energies  and  vrave  vectors  that  satisfy 

|E  - Wa(k„)  + eAv|  -2Ta(k„) 

|E  - Wb(k„)  - eAvl  >2Tb(k„) 
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the  corresponding  wave  functions  are  Bloch-like  inside  metal  a but  decay  exponen- 
tially from  the  interface  inside  metal  b . In  a similar  way,  the  wave  functions 
of  states  whose  energies  and  wave  vectors  satisfy 


iE  ' wa  Ql.i)  + eAv|  >2T  (k„) 


|E  - Wb(k„)  - eAv|<2Ta(k„) 


(4.3) 


extend  to  infinity  on  the  b side  of  the  interface,  but  decay  exponentially 
on  the  a side.  The  decay  coefficient  of  the  wave  function  is  determined  by 
the  energy  measured  relative  to  the  corresponding  k,,  sub-band  center.  If  we 
express  this  energy  in  units  of  the  sub-band  half -width,  and  write 


(k„)  = o2T(k„)  , M 1 


(4.4) 


where  to(k  ,,)  and  T(k,,)  refer  to  the  corresponding  values  on  the  side  of  the 
interface  where  the  wave  function  decays  exponentially,  we  obtain  from  Eqs. 
(3.4)  and  (3.23)  that  on  the  decaying  side  the  wave  function  on  plane  m is 
proportionally  to  e '>Jn  . The  decay  coefficient  X>0  is  given  by 


X - - In  I 


sign (a) 


y *2-i] 


(4.5) 


Thus,  the  further  the  energy  is  from  the  center  of  the  sub-band,  the  stronger 
is  the  exponential  decay  of  the  corresponding  wave  function. 

We  note  that  the  behav.^,.  of  the  wave  functions  extending  througliout  the 
entire  system  and  those  extending  only  on  one  side  of  the  interface  is  as  expected. 
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If  we  try  to  propagate  a wave  through  the  interfaced  crystal  this  wave  can  prot_- 
gate  from  one  side  to  the  other  only  it  its  frequency  is  in  the  cannon  sub-bands. 
If,  however,  the  frequency  is  in  the  sub-band  of  one  of  the  crystals  but  outside 
the  sub-band  of  the  other,  this  wave  cannot  penetrate  into  the  second  crystal, 
and  its  amplitude  will  decay  exponentially. 

The  third  type  of  wave  function  is  associated  with  the  existence  of  bound 
interface  states.  For  a given  k,,,  the  energy  of  the  possible  bound  states  are 
determined  by  the  poles  of  the  Green's  function,  which  lie  outside  the  k,,  sub- 
bands of  the  two  metals  forming  the  junction.  Using  the  explicit  expression 
(3.23)  for  the  interface  Green's  function,  we  see  that  the  bound  states  energies 
are  given  by  the  roots  of  the  equation: 

402  - ^sign(wa)fla+  wa  - 2u|ign(u>b)jjb+  ^ - 2uJ  = 0 5 

outside  the  k,  sub-bands.  In  this  equation  y is  defined  by 


’ = (u)2  - 4T2)^ 

As  can  be  seen  from  Eq.  (3.23),  the  wave  function  of  an  interface  bound  state 
is  localized  near  the  interface  and  decays  exponentially  with  distance  on  both 
sides  of  the  interface.  Thus,  the  bound  electron  is  free  to  propagate  parallel 
to  the  interface,  but  is  confined  to  a finite  stripe  of  width  A^  +A^  (where 
A and  A,  are  given  by  Eq . (4 . 9)  below)  in  the  direction  perpendicular  to  the 
interface.  The  decay  coefficients  of  the  bound  state  wave  function  will  be 
different,  in  general,  on  the  two  sides  of  the  interface.  If  we  express  the 


bound  state  energy,  relative  to  the  center  of  the  two  k,.  sub-bands,  in  terms  of 
the  corresponding  sub-bands  half-widths , we  can  write 


< = aa2Ta,  l*J>i  ; < = *b  2V  KI>1 


(4.7) 


where 


U)°  = E°  - W + eAv  : = E°  - W.  - eAv 

a a ’ b b 


(4.3) 


and  E°  is  the  energy  of  the  bound  state.  It  is  easy  to  see  frcm  Eq.  (3.23)  that 
the  exponential  decay  coefficients  on  side  a and  side  b , and  A^  res- 


pectively, are  given  by 


a = - ln[aa  - Ja a2'1  ] 

‘n  = “ in[»b  ~ signCQ^^V^2-1  ] 


[®a  ' s ign  (®a}  'JaB 


(4.9a) 


(4.9b) 


Thus,  the  greater  the  distance  of  the  bound  state  energy  frcm  the  center  of  the 
k,,  sub-bands,  the  more  localized  is  the  corresponding  wave  function. 

We  note  that  a general  property  of  the  bound  interface  states,  derived 
from  Eq.  (4.6),  is  that  whenever  two  interface  states  exist  simultaneously,  an 
increase  in  the  coupling  constant  6 will  increase  the  energy  of  the  interface 
state  having  the  higher  energy  and  will  decrease  the  energy  of  the  lower  energy 
bound  state.  Another  remark  that  should  be  added  here  is  that  for  very  strong 
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coupling  between  the  two  metals  there  will  always  be  two  interface  states,  a bonding 
state  be?ow  the  sub-bands  and  an  anti-bonding  one  above  the  sub-bands.  The 
asymptotic  energies  of  these  bound  states  are  given  by 

E1  = Ua  + wa  " eAv  + 3 
E2  = ^b  + ^b  + eAv  ” ^ 

if  j Ua  + - eAv|  > |U^  4-  + eAv  | and  by 

E-^  — + eAv  + 6 

E0  = U + W -eAv  - 6 

Z a a 

if  I + ^b  + e^v  I > | + Wa  - eAv  | 

The  analytical  solution  of  Eq.  (4.6)  for  the  bound  states  energies  is  not  possible 
in  general.  However,  by  a suitable  graphical  analysis,  one  can  determine  the 
conditions  for  the  occurance  of  bound  states,  their  number  and  their  position 
with  respect  to  the  bands,  as  a function  of  the  strength  of  the  coupling  constant. 

It  turns  out  to  be  very  convenient  to  analyze  the  interface  bound  states  in 
terms  of  the  properties  of  the  bound  surface  states  of  the  intermediate  surface 
problem  considered  in  the  preceeding  section  (i.e.  the  surface  problem  having  a 
surface  self-consistent  potential  equal  to  the  respective  value  of  the  interface 


7‘r  ' - 


3 


problem) . It  is  well  known  that  the  bound  surface  states  are  given  by  the  roots 
.(13). 


of  the  equation^ 

sign(co)  y -+w  -2U=  0 

The  corresponding  surface  bound  state  energies  are  given  by 


o Ta  (— 

E - W (k„)  - eAv  + U + — ~ 

a a — a Ua 


(4.10a) 


Eb  - Wb(k„)  + eAv  + Ub  + 


Tb2  <&,) 


(4.10b) 


provided  |Ua|>Ta  and  | | >T^  . For  U < T there  are  no  surface  states. 

In  the  following  analysis  we  shall  use  the  indices  1 and  2 to  denote  the  metal 
having  the  higher  sub-band  top  edge  and  the  lower  sub-band  bottom  edge,  respec- 
tively. 

The  behavior  of  the  interface  bound  states  can  be  described  as  follows: 
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and 


fi(E)=sign(0)i)tIi+  aii-2Ui 


(A. 12) 


E-^,  max  and  E2  min  are  the  top  edge  and  the  bottom  edge  energies  of  the  higher 

and  the  lower  of  the  two  k,,  sub-bands  respectively.  If  -(T2-HJ2)f^(E2  j^)  < 

(T,  - Ui)f0(E1  ) the  bound  state  will  appear  below  the  sub-bands.  If.  how- 

1 1 / l,max 

ever,  -(T2+U2)f1(E2  uax)  > (T^  - U1)f2(E1  max)  the  bound  state  will  appear  above 
the  sub- bands.  The  bound  state  described  above  exists  as  long  as  8>B^  . When 
the  coupling  constand  is  further  increased  and  the  value  B2  is  reached,  a 
second  bound  state  appears.  The  second  critical  coupling  is  given  by 


,1 


max 


[- 


%(T2+U2)fl(E 


).  %(T  -U  )f  (E 
2 min  I .max' 


1 \ 

J( 


(4.13) 


for  $ > £>2  there  are  always  two  bound  states,  one  below  and  the  other  above 
the  sub-bands. 


b.  Suppose  that  in  the  intermediate  surface  problem  one  of  the  following 
four  situations  holds: 

1.  There  is  one  bound  state  above  the  sub-bands  and  the  two  sub-bands 
overlap . 

2.  There  is  one  bound  state  above  the  sub-bands,  which  is  associated 
with  the  metal  having  the  higher  sub-band  and  there  is  a gap  between  the 
two  sub-bands. 
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3.  There  are  two  surface  states,  one  above  the  sub-bands  and  the  other 
lies  in  the  sub-band  of  the  second  metal  and  there  is  an  overlap  between 
the  sub-bands. 

4.  There  are  two  bound  states,  a gap  exists  between  the  sub-bands  and 
the  bound  state  of  the  lower  sub-band  lies  above  the  sub-bands,  whereas, 
the  one  associated  with  the  upper  sub-band  falls  in  the  lower  sub-band. 


Then  the  behavior  of  the  interface  bound  states  is  as  follows:  for  small  coupling 
constant  there  will  be  one  bound  state  above  the  sub-bands.  This  state  develops 
in  a continuous  way  from  the  corresponding  surface  state.  If  E?  is  the  energy 
of  the  surface  state,  which  lies  above  the  sub-bands,  then  for  small  g the 
energy  of  the  corresponding  interface  state  is  given  approximately  by 

32 

_ n - Si  > 

(4.14) 


E ~ E° 

— l 


O d - T > 

f J 


where  f^  (E)  is  the  function  defined  by  (4.12).  As  long  as  6 is  in  the  range 


o<£<62  where  £$3  is  defined  by 

= , = -%  (T,+U„)f  e . ) 
3 {_  ill  2, mm 


(4.15) 


there  is  one  bound  state  above  the  sub-bands.  For  g^g^  there  are  two  bound 
interface  states,  one  above  and  the  other  below  the  sub-bands. 

c.  Suppose  that  in  the  intermediate  surface  problem  one  of  the  following 
four  situations  exists: 
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1.  There  is  an  overlap  between  the  sub-bands,  and  a bound  state  exists 
below  the  sub-bands. 

2.  There  is  a gap  between  the  sub-bands,  there  is  one  surface  state 
below  the  sub-bands  which  is  associated  with  the  lower  sub-band. 

3.  There  is  an  overlap  between  the  sub-bands  and  there  are  two  sur- 
face states,  one  below  the  sub-bands  and  the  other  lies  inside  the 
sub-band  of  the  second  metal. 

4.  There  is  a gap  between  the  sub-bands  and  there  are  two  surface 

states,  the  one  associated  with  the  higher  sub-band  lies  below  the  sub-bands, 
and  the  other  falls  inside  the  higher  sub-band. 


Then  the  behavior  of  the  interface  bound  states  is  as  follows:  for  small  6 there 
will  be  one  bound  state  below  the  sub-bands,  with  an  energy  given  approximately 
by  (4.14).  When  the  coupling  constant  is  increased,  a critical  value  6^,  is 
reached,  where  a second  Lxjund  state  appears  above  the  sub-bands,  g^  is  defined  by 


(4.16) 


For  g>_g^  there  are,  thus,  two  bound  states,  one  below  and  the  other  above  the 


sub-bands . 


d.  Suppose  that  in  the  intermediate  surface  problem  one  of  the  following 
situations  holds: 

1.  There  is  a gap  between  the  sub-bands  and  there  exists  one  surface 
state  above  the  sub-bands,  which  is  associated  with  the  lower  sub-band. 
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2.  There  is  a gap  between  the  sub-bands  and  there  exist  tvro  surface 
states.  The  one  associated  with  the  higher  sub-band  is  located  above 
the  sub-bands  and  the  other  lies  inside  the  higher  sub-band. 


Then  the  behavior  of  the  interface  states  is  as  follows:  for  a small  coupling 
constant  there  is  one  interface  state  above  the  sub-bands.  When  8 reaches  the 


critical  value  3 where 
5 


-U  tt1  + 


> f2  <E,  . ) 
l,min. 


(4.17) 


(E,  . being  the  bottom  edge  energy  of  the  higher  sub-band)  a second  bound 

v l,min 

state  appears  inside  the  gap.  This  bound  state  exists  as  long  as  the  coupling 
constant  is  in  the  range:  8^f_  8 £ 8^,  where  8g  is  given  by 


— I 


h (T0  - U„)  f.  (E 


) 


(4.18) 


2, max 


E„  is  the  energy  at  the  top  of  the  lower  sub-band.  For  8>8  r the  bound 
2,max  ° 

state  disappears  from  the  gap.  When  8 reaches  the  value  given  by  Eq. 

(4.15),  a new  interface  state  appears  below  the  sub-bands.  For  g>g3  there 
are  two  interface  states,  one  below  and  the  other  above  the  sub-bands. 

e.  Suppose  that  in  the  intermediate  surface  problem  one  of  the  following 
situations  holds: 
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1.  There  is  a gap  between  the  sub-bands,  and  there  exists  one  surfac- 
state  associated  with  the  higher  sub-band  which  lies  below  the  sub-bands. 

2.  There  is  a gap  between  the  sub-bands  and  there  are  two  surface  states 
The  one  associated  with  the  lower  sub-band  lies  below  the  sub-bands  and 
the  other  one  lies  in  the  lower  sub-band. 

Then  the  behavior  of  the  interface  states  is  as  follows:  for  small  values  of  the 
coupling  constant  there  is  one  interface  state  below  the  sub-bands,  whose 
energy  is  given  approximately  by  (4.14).  When  6 is  increased,  and  reaches 
the  value  given  by  (4.18),  a new  interface  state  appears  in  the  gap.  This 
bound  state  exists  in  the  gap  for  B^<B  £85  where  B^  is  given  by  (4.17). 

For  S > 65  the  bound  state  disappears  from  the  gap  and  a new  bound  interface 
state  appears  above  the  sub-bands  for  3 £ 3^>  where  8^  is  given  by  (4.16) . 

For  6 > 84  there  are  two  bound  states,  one  below  and  the  other  above  the  sub- 
bands. 

f.  Suppose  that  one  of  the  following  situations  holds  in  the  intermediate 
surface  problem: 

1.  There  is  a gap  between  the  sub-bands  and  there  is  a single  surface 
state  which  lies  in  the  gap. 

2.  There  is  a gap  between  the  sub-bands  and  there  are  two  surface 
states,  one  in  the  gap  and  the  other  lies  inside  the  sub-band  of  the 
other  metal. 

Then  the  behavior  of  the  interface  bound  state  is  as  follows : for  small  values 
of  the  coupling  constant  there  is  a single  interface  state  which  lies  in  the 


gap.  When  g is  increased  and  readies  the  value  By  where 


n \ % 


(4.19) 


By  =|max  [^(T^)  f2  (E^^),  % (VV  fj_  (E2inax)  j 

this  bound  state  disappears.  For  8>8  there  appears  a new  interface  state 

3 

below  the  sub-bands  where  83  is  given  by  (4.15).  If  8>B4  there  is  also 
a bound  state  above  the  sub-bands , where  B4  is  determined  by  (4 . 16) . 


g.  Suppose  that  in  the  intermediate  surface  problem  there  are  two  sur- 
face states  below  the  sub-bands,  and  a gap  exists  between  these  sub-bands, 
then  for  small  values  of  B there  will  be  two  interface  states  below  the 
sub-bands  (with  energies  given  approximately  by  (4 . 14) ) . When  8 reaches 

the  value  given  by  (4.15),  the  interface  state  having  the  higher  energy  disap- 
pears . A new  bound  interface  state  appears  in  the  gap  when  B reaches 

value  of  8,  where  BA  is  given  by  (4.18).  This  bound  state  exists  as 
6 o 

lo»ig  as  B^B£B5  where  85  is  determined  by  (4.17).  For  B > 85  the 
interface  state  disappears  from  the  gap  and  a new  bound  state  appears  above 
the  sub-bands  when  B _>  B4  , where  84  is  given  by  (4 . 16) . 

h.  Suppose  that  the  situation  in  the  intermediate  surface  problem  is  the 
same  as  described  in  (g.)  except  that  there  is  an  overlap  between  the 
sub-bands.  Ihere  will  be  two  interface  states  below  the  sub-bands  for 

0 < B < 83,  where  83  is  given  by  (4.15).  A new  interface  state 
appears  above  the  sub-bands  for  6 > 8^  , where  B4  is  given  by  (4 . 16) . 
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i.  Suppose  tliat  in  the  intermediate  surface  problem  there  are  tvro  surface 
states  above  the  sub-bands,  and  a gap  exists,  then  for  srv-11  values  of  3 
there  are  two  interface  states  above  the  sub-bands.  When  8 ’reaches  the  value  8^ 
given  by  (4.16),  the  interface  state  having  the  lower  energy  rusappears.  A 
new  bound  state  appears  inside  the  gap  when  3 reaches  the  varre  8^  given  by 
(4.17).  This  bound  interface  state  exists  as  long  as  85  £ 8 i % where 

g is  given  by  (4.18).  If  3 is  further  increased,  the  bounc  state  disappears 
from  the  gap  and  a new  interface  state  appears  below  the  sub-hinds  when  B 
where  83  is  given  by  (4 . 15) . 

k.  Suppose  that  in  the  intermediate  surface  problem  tnere  are  two  surface  states 
inside  the  gap,  then  for  snail  values  of  the  coupling  constant  there  will 
be  two  interface  states  inside  the  gap.  Their  energies  are  given  approxi- 
mately by  (4 . 14) . When  8 reaches  the  value  8g  where 


“8  -j"1"  b <VV  fl<E2„«).  -4  <T1  + V r2  j 


(4.20) 


one  of  these  bound  states  disappears.  When  8 is  further  increased  and 
reaches  the  value  89  given  by 

| % 

N -U[  ‘l  (E2iJ'  ^ <T‘  + ^ "i.Jh 


(4.21) 
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The  other  bound  states  also  disappears  from  the  gap.  For  B ■>  B3.  where 
B3  is  given  by  (4.15),  an  interface  state  appears  below  the  sub-bands, 
whereas  for  B^  B^  , where  B^  is  given  by  (4 . 16) , there  is  a bound  state 
above  the  sub-bands. 

l.  Suppose  that  in  the  intermediate  surface  problem  there  are  two  surface 
states,  one  below  and  the  other  above  the  sub-bands,  then  for  any  value  of 
the  coupling  constant  B there  will  be  two  interface  states,  one  below  and 
the  other  above  the  sub-bands. 

m.  Suppose  that  in  the  intermediate  surface  problem  there  are  two  surface 
states,  one  in  the  gap  and  the  other  below  (above)  the  sub-bands,  then  for 
small  values  of  B there  will  also  be  two  interface  states,  one  in  the  gap 
and  the  other  below  (above)  the  sub-bands.  The  bound  state  in  the  gap 
disappears  for  B > By  , where  By  is  given  by  (4 . 10) . For  6 _>  B^ 

(B  i 63)  a new  bound  state  appears  above  (below)  the  sub-bands,  where 
B^  (B3)  is  given  by  (4.16),  and  (4.15)  ). 

This  analysis  covers  the  various  possible  bound  interface  states. 

In  Section  VI  we  shall  return  to  the  problem  of  the  interface  states  while 
discussing  a numerical  example. 

V.  The  Interface  Energy 

If  many-body  effects  are  neglected  the  total  energy  of  the  crystal  with 
the  interface  is  just  the  sum  of  the  occupied  single-particle  energy  levels. 
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This,  in  turn,  can  be  expressed  as  a corresponding  integral  over  the  system's 
density  of  states.  The  total  density  of  states  of  the  interfaced  crystal  can  be 
obtained  from  the  Green's  function  derived  in  Section  III,  by  sunming  the 
imaginery  part  over  the  various  planes  parallel  to  the  interface.  However,  in 
order  to  determine  the  energy  needed  to  break  the  metal -metal  interface  into  two 
metal-vacuum  interfaces  we  need  to  know  the  change  in  the  total  density  of  states 
due  to  this  cleavage  process.  This  can  be  evaluated  directly,  without  calculating 
the  density  of  states  of  the  two  systans,  by  using  the  phase  shifts  method  des- 
cribed by  De  Witt^^,  Callaway  and  Toulouse  . This  method  can  be 

surrmarized  as  follows:  let  H = H°  + V be  the  perturbed  Hamiltonian,  and  let 
and  be  the  eigenvalues  of  H and  H°  respectively.  The  change  in  density 
of  states  due  to  the  perturbation  V can  be  written  as 


Ap  (E)  = p (E)  - o°  (E)  = E j^5  (E-E^ ) - (E-eu)] 


where  p and  p are  the  unperturbed  and  the  perturbed  densities  of  states. 
This  can  be  also  expressed  in  the  following  way: 


Ap(E)  = 


T , r E-i5-E . 

_ Im  o 1 

' ~ Se  z . . . 

j E - l C- o . 

*-  2 

■ h ST  Tr  °'lc°] 


where  G and  G°  are  the  Green's  operators  of  the  perturbed  and  the  unperturbed 
Hamiltonians,  respectively.  Using  Dyson's  equation  we  can  express  G-^  in 
terms  of  Gr  . In  this  way  we  get 
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Ao(E)  Im^Y  Tr  j^n  (1-G°V)J 
= n Imaf  in  [det  <1-G°V)] 


(5.1) 


In  the  last  step  we  used  the  identity  Tr  JfnA  J = fn  Jdet  A j , where  A is 
a general  operator.  The  application  of  Eq.  (5.1)  is  especially  useful  for  a 
localized  perturbation,  when  det  (1-G°V)  can  be  expressed  as  a finite-order 
determinant . 

Let  us  turn  now  to  the  specific  problem  of  the  interface.  In  this  case  we 
start  from  the  two  semi-infinite  surface  systems  (the  a and  the  b surfaces), 
and  apply  the  perturbation  necessary  to  create  the  interface.  This  perturba- 
tion was  described  in  detail  in  Section  V.  We  introduce  the  perturbation  in  two 
steps.  In  the  first  one  we  apply  a constant  electrostatic  potential  Av  on 
metal  a and  -Av  on  metal  b . As  we  have  seen  before,  the  application  of 
these  potentials  causes  a shift  in  the  densities  of  states  of  the  two  surface 
systems , which  aligns  the  Fermi  levels  of  the  two  metals.  Since  each  metal  is 
electrically  neutral,  there  will  be  no  net  change  in  their  energies  due  to  the 
application  of  this  perturbation.  Our  second  and  final  step  is  to  apply  the 
remaining  perturbation  V needed  to  form  the  interface.  As  we  have  seen  in 
Sections  II  and  III,  the  only  non- vanishing  matrix  elements  of  this  perturba- 
tion are  given  by 


V (0,0)  = ua  - ua 

(5.2a) 

V (-1,-1)  ub  - Ub 

(5.2b) 

» 


-33- 


......  ■ . -L  . > " ~ ^ 


and 


v (-1,0)  = V(0,-1)*  = 


i ^ 


(5.2c) 


where,  as  before,  Ua>  U^,  and  U^,  0^  are  the  self-consistent  potentials  near 
the  interface  and  near  the  surface,  respectively.  Thus,  in  the  present 
case  we  have 


det  (1-G  V)  = 


1 - G (-1 , -1)  V (-1,-1) 


-G  (0,0)  V (0-1) 


-G°  (-1,-1)  V (-1,0) 


1 - G (0,0)  V (0,0) 


(5.3) 


where  G°(0,0)  and  G°(-l,-l)  are  the  corresponding  surface  Green's  functions 
of  metals  a and  b,  respectively,  with  the  energies  shifted  according  to  (3.17). 
We  note  that  the  surface  self-consistent  potentials  in  G°(0,0)  and  G°(-l,-l) 
are  and  0^  , respectively,  and  not  and  , as  was  the  case  in  the  pre- 
ceding section.  Since  G°  and  V are  diagonal  in  the  wave  vector  k,,  , the  deter- 
minant (5.3)  factorizes  into  similar  terms  with  different  k,,  values.  Expanding 
(5.3)  in  terms  of  the  indices  m and  n shows  that  each  such  factor  is  given 


by 


detk„  (1-G°V>  = [l-v<-l,-D  G°(-l,-l)][l-V(0,0)  G°(0,0)j  -e2G°  (0,0)G°(-1,-1) 


(5.4) 


In  this  expression  we  used  the  notation  det^  for  the  partial  determinant,  with 
a specific  k,,  . 
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Substituting  the  explicit  expressions  for  the  surface  Green's  functions  G°(0,0) 
and  G°(-l,-l),  from  the  general  expression  (3.21),  into  Eq.  (5.4)  gives 

_ 4P2-fk+1(Xa-2Ua)1  k^v21^ 


det,  (1-G  V) 


(5.5) 


The  bound  state  energies  of  the  crystal  with  the  interface  are  given  by  the 
roots  of  the  equation^  ^ : 


det  (1-G  V)  = 0 


(5.6) 


outside  the  shifted  bands  of  the  two  metals.  Using  the  explicit  expression  (5.5) 
we  see  that  this  is  exactly  the  same  condition  derived  earlier  from  the  poles 
of  the  interface  Green's  function,  Eq.  (4.6). 


shifts  by 


In  analogy  with  ordinary  scattering  theory,  one  defines  the  partial  phase 
(16-17) . 


T]  (E,k„)  = Tm  In 

!detk„  (l-C^M 

L *“  J 

= Arg  | 

|4B2+  L + i<v2V 

k + i<v2v 

1 1 

a + i(“>  - 2 0 )1 

a a a 

-b  + 1 <V2V 

(5.7) 


In  the  second  step,  we  applied  the  explicit  expression  (5.5)  and  used  the  iden- 
tity Imtnf  = Argf , where  Arg  denotes  the  argument  of  the  complex  function. 

If  we  take  the  determinant  of  (5.5)  with  respect  to  k,.,  and  use  the  rela- 
tion (5.1)  we  see  that  the  change  in  the  density  of  states  per  surface  atom, 
relative  to  the  shifted  bands,  due  to  the  formation  of  the  interface,  is  given  by: 
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-.  tWtfdLafc, 


A P<E)  • $ 72^)2  If  / <%"  H®*’ 


all- 

TT  (2TT)'1 


( 4B2+  [ua  + i(“a  " 2tVl  MWlj, 

_ /dk„  Arg  j I + , («,a-2Ua)f(II7T  1 <V2fib>l  ' 


In  this  expression  we  introduced  explicitly  the  g-fold  degeneracy  of  the  bands 
under  consideration.  We  also  note  that  An  is  the  area  of  the  unit  cell  para- 
llel to  the  interface  and  that  the  k,.  integration  goes  over  the  two-dimen- 
sional Brillouin  zone  defined  by  the  crystals  structure  parallel  to  the  interface 
Equation  (5.8)  can  be  used  to  derive  a relatively  simple  expression  for  the 
change  in  the  integrated  density  of  states.  This  is  given  by: 


AN  (E)  = 


Ap(E)  dE 


{ 4B2+  i (V2V  + l(V2V) 

- / dk„  Arg<-f 7-T ~ ] ( 

)2  J \ (j,  + i(w  -2U  ' - i(%-2Ub) 


The  single-particle  contribution  to  the  interface  energy  at  T=0  (i.e. 
the  energy  per  interface  atom  necessary  to  break  the  interface  into  tws  semi- 
infinite  crystals,  at  T=0)  is  given  by 
E’ 

iab  = - J A p(E)  EdE  (5-1 

“oo 

where  the  Fermi  energy  Ej  differs  from  the  value  given  by  Eq.  (2.2)  by  a term 
which  is  0(1/Nj) . This  E£  guarantees  the  charge  neutrality  of  the  interfaced 


- 


crystal.  Expanding  (5.10)  to  first  order  in  Ej  - Ff  we  can  express  the  inter- 
face energy  as  follows: 

E* 

(5.11) 


ab 


/ 


o.v  = - / Ap(E)(E-Ef)  dE 


■ L 


AN(E)  dE 


, £.  Jl- 

" (2TT)2 


4c>2+ 

K.  + i(w, 

a a 

- 2U  ) 
a 

M,b+L(V 

• + 
a 

i(u)  -2U  ) 
a a 

%+i(a,h2eb) 

sl  fo 

Also,  is  the  common  Fermi  energy  of  the  tro  metals,  i.e.  E^  = % (E^  + E ) and 

ooa  and  include  the  corresponding  shifts  in  the  energies,  according  to  Eq. 
(3.22). 

In  a similar  way,  one  can  determine  the  single-particle  contribution  to  the 
change  in  the  electronic  specific  heat.  This  change  is  given  by: 

cc 

ACV  = / EAp(E)  ||  dE 

JC 

where  f(E)  is  the  Fermi-Dirac  distribution  function.  For  temperatures  much  lower 
than  the  Fermi  temperature  the  change  in  the  electronic  specific  heat  is  linear 
in  the  temperature  and  proportional  to  the  change  in  the  density  of  states  at 
the  Fermi  level  (assuming  no  Van  Hove  singularity  occurring  at  the  Fermi  energy) . 
Explicitly  we  have 
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where  Kg  is  the  Boltzmann  constant.  For  the  interface  system  the  constant  y is 
given  by 


L2n 


A„ 


/d-"  jE 


Arg 


+ 

-2  + i(V 

-2U  ) 
a 

M,  +1(01 

D b 

l°v] 

pu  + 

a 

i(u>  -2U  ) 
a a 

> 

1 <V2V 

(5.12) 


where  the  integrand  has  to  be  evaluated  at  the  comnon  Fermi  level. 

VI.  Numerical  Results  and  Discussion 

As  an  application  of  the  formalism  developed  in  the  last  three  sections  we 
consider  in  the  following  the  interface  formed  parallel  to  the  (100)  plane  of  two 
crystals  described  by  the  same  tight-binding  parameters,  but  having  different 
Fermi  energies.  The  coupling  between  the  two  metals  is  taken  to  be  the  same  as 
the  bulk  coupling,  we  also  neglect  the  interface  perturbation  and  (i.e. 
we  sec  Ta  = Tg  = B , and  = 0) . This  simplification  is  done  here  in 

order  to  reduce  the  number  of  independent  paraments  in  the  problem.  The  general 
case,  discussed  earlier,  can  be  analyzed  in  a completely  analoguous  way.  The 
specific  case,  considered  here,  can  serve  as  a crude  model,  describing  the  elec- 
tronic properties  of  a (100)  interface  formed  between  two  transition  metals 
belonging  to  the  same  series,  and  will  be  referred  to  as  such  in  the  followring. 
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From  the  preceeding  discussion,  presented  in  Section  IV,  it  is  obvious 
that  there  are  no  interface  bound  states  in  the  model  under  consideration.  The 
electronic  wave  functions  are  delocalized,  and  extend  on  either  one  side  or  on 
both  sides  of  the  interface,  according  to  the  corresponding  electron  energy.  For 
the  (100)  interface  we  have 


(k,,)  = 2E^  J^cos(ak + cos(akz)j 


T = E, 


(6.1a) 


(6.1b) 


Using  these  relations  and  Eq.  (3.23)  for  the  diagonal  matrix  elements  of  the 
Green's  function,  and  the  expression  (3.11),  we  can  calculate  the  local  densities 
of  states  on  the  various  planes  of  the  interfaced  crystal.  Figure  (1)  shows  the 
results  of  such  a numerical  calculation  of  the  local  density  of  states  pn(E) , for 
the  first  three  atomic  layers  adjacent  to  the  interface  (n  = 0,  1 and  2,  respec- 
tively) . The  difference  in  the  Fermi  energies  of  the  two  metals  was  chosen  to 

be  AE^  = 2E^  . Fig.  (2)  shows  the  corresponding  density  of  states  for  the  case 

of  AE^  = 8E^.  For  comparison,  the  shifted  bulk  density  of  states  is  also  shown 
in  these  figures.  All  the  density  of  states  presented  here  are  normalized  to 
unity  (i.e.  we  set  g = 1) . We  note  that  Figs.  (1)  - (2)  refer  to  the  side  of 
the  interface  on  which  the  metal  having  the  higher  Fermi  energy  is  located.  The 
local  densities  of  states  on  the  other  side  of  the  interface  can  be  obtained  from 
the  curves  in  Figs.  (1)  and  (2)  by  taking  their  mirror  images,  with  respect  to 

the  E = 0 line.  As  can  be  seen  from  the  above-mentioned  figures,  the  local  densi- 

ties of  states  have  a "tail"  extending  outside  the  shifted  bulk  band.  This  tail 
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is  contributed  entirely  by  electrons  tunneling  from  the  other  side  of  the  int'  -face. 
From  Figs.  (1)  and  (2)  it  is  obvious  that  the  penetration  distance  of  these  elec- 
trons is  essentially  limited  to  only  a few  atomic  layers.  In  general,  the  trost 
pronounced  effect  on  the  local  density  of  state  occurs  near  the  interface  itself 
(i.e.  n = 0).  As  one  proceeds  away  from  the  interface,  the  density  of  states 
approaches  asymptotically  the  bulk  density  of  states.  For  n = 2 (the  third 
layer  near  the  interface)  the  density  of  states  of  the  interfaced  crystal  is 
already  very  close  to  the  corresponding  shifted  bulk  density  of  states. 

When,  in  the  model  under  consideration,  the  coupling  constant  between  the 
two  metals  is  allowed  to  vary,  it  is  possible  to  form  bound  interface  states. 
Following  the  discussion  of  Section  IV  it  is  not  difficult  to  show  that  for 
the  present  model  the  two  critical  values  B^  and  coincide.  Therefore, 
in  this  case  there  will  be  either  no  bound  states  at  all  or  there  will  be  two 
of  them,  one  above  and  the  other  below  the  sub-bands.  It  can  be  easily  shown  that 
the  critical  coupling  constant  Bc,  to  which  both  and  reduce,  is  related 

to  the  difference  in  the  Fermi  energies  of  the  two  metals  as  follows 


6C  = I T + % AEf  +*/%  AEf  + T AEf 


(6.2) 


Thus,  the  bigger  the  difference  in  the  Fermi  energies  the  larger  is  the  coupling 
necessary  to  create  a bound  state.  This  is  quite  expected  since  the  greater  is 
the  difference  in  the  Fermi  energies  the  larger  is  the  minimum  distance  of  the 
bound  states  from  the  center  of  the  sub-bands.  The  minimum  possible  critical 
coupling  is  obtained  for  AE^  = 0 , and  is  given  by  Bc  = T.  Therefore,  in 
the  present  model,  every  coupling  constant  which  exceeds  the  bulk  coupling( vri.ll 
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produce  a bound  state.  Fig  (3)  shows  the  dependence  of  the  bound  state  energies 
on  the  coupling  constant  for  the  case  of  AE^  = 0 and  AE^  = 2T.  As  can  be 
from  these  curves,  an  increase  in  the  coupling  constant  increases  the  energy 
of  the  upper  bound  state  and  decreases  that  of  the  lower.  This  behavior  agrees 
with  our  general  discussion  given  in  Section  IV. 

We  now  turn  to  the  calculation  of  the  change  in  the  total  density  of  states 
and  the  interface  energy  of  two  transition  metals  belonging  to  the  same  series. 
Using  Eq.  (5.8)  it  is  not  difficult  to  see  that  the  change  in  the  density  of 
states,  in  the  model  under  consideration,  is  an  even  function  of  the  energy  rela- 
tive to  the  center  of  the  band.  Using  Eq.  (5.9)  we  calculated  numerically  the 
change  in  the  integrated  density  of  states  AN(E)  due  to  the  creation  of  the 
interface.  Fig.  (4)  shows  the  change  in  the  integrated  density  of  states, 
which  is  an  odd  function  of  the  energy,  for  the  cases  of  AEf  = 2E^  , AEf  = 6E^ 
and  AE^  = 12E^,  respectively.  In  this  figure  we  used  the  value  g=10  to  account 
for  the  10- fold  degeneracy  of  the  d orbitals  in  transition  metals.  To  obtain 
the  corresponding  interface  energy  we  have , according  to  Eq . (5 . 11) , to  integrate 
AN(E)  with  respect  to  the  energy,  up  to  the  cannon  Fermi  energy.  The  results 
of  such  a numerical  calculation  are  shown  in  Fig.  (5),  where  the  interface  energy 
for  a given  transition  metal  with  n^  d electrons  is  ploted  vs  , the  number 
of  d electrons  of  the  other  metal  forming  the  interface.  The  curves  in  this 


figure  correspond  to  the  cases  where  nfl  = 0,  1,  . . .5.  The  corresponding  graphs 
for  ^>5  (i.e.  = 6,  ...10)  can  be  obtained  by  taking  the  mirror  images  of 
the  curves  of  10  - na  electrons,  with  respect  to  the  = 5 line.  This  is 


due  to  the  face  that  in  the  present  model,  the  following  symnetry  holds 


-41- 


a(na,  r^)  = a(10  - 10  - r^) 

where  <5  and  n^<5.  The  surface  energy  of  transition  metals  can  be  read 
off  the  curves  of  Fig.  (5)  by  looking  at  the  points  where  • Our 

results  for  the  surface  energies  are  the  same  as  those  derived  earlier  by  F. 

Cyro t -Lackmann and  by  G.  Allan.  This  is  as  expected,  since  in  the 

limit  of  two  identical  transition  metals  our  interface  model , discussed  in  this 
section  reduces  to  the  free  surface  model  investigated  by  these  authors . 

As  can  be  seen  from  the  Fig.  (5),  the  electronic  contribution  to  the  inter- 
face energy  between  two  transition  metals,  of  the  same  series,  is  always  less 
than  the  corresponding  contribution  to  the  surface  energy  of  a half -filled 
band  metal  (belonging  to  the  same  series) . We  also  note  that  for  a given  transi- 
tion metal,  with  d-electrons,  there  exists  another  transition  metal,  with 
n^°  electrons,  whose  combined  interface  has  the  maximum  interface  energy.  As 
na  is  varied,  the  position  of  this  maximum  drops  down  from  n^°~8  for  an  empty 
band  (na  = 0)  to  n^°  ~2  for  a full  band  (n^  = 10) . 

The  interface  model  developed  in  the  present  work  is  a highly  simplified 
and  crude  one.  Nevertheless,  we  believe  that  certain  features  of  a real  inter- 
face, such  as  interface  states  and  electron  tunneling  across  the  interface, 
are  illustrated  by  our  model.  Thus,  although  a realistic  physical  description 
of  a bimetallic  interface  might  be  much  more  complicated,  we  feel  that  many  of 
the  qualitative  properties  of  such  an  interface  will  be  similar  to  those  des- 
cribed in  this  article. 
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APPENDIX 

In  this  Appendix  we  derive  the  expressions  for  the  general  matrix  elements 

of  the  surface  Green's  function.  The  technique  we  apply  is  the  same  as  that 

(13) 

used  by  Kalkstein  and  Soven  to  derive  the  corresponding  diagonal  matrix 
elements . 

To  form  to  metal-vacuum  interface  we  start  from  the  bulk  metal  and  break 
the  bonds  between  the  tvro  adjacent  planes  n = 0 and  n = -1.  In  addition, 
we  change  the  self-consistent  potential  of  the  electrons  near  the  surface.  The 
perturbation  described  above,  which  is  necessary  to  create  the  surface,  has 
the  following  non-vanishing  matrix  elements  in  the  Bloch-Warmier  representation: 


v(0,0)  = v (-1,-1)  = u 


V(-1,0)  = V(0,-1)*  = -Te 


10 


(A.  1) 


(A.  2) 


With  this  perturbation  potential,  the  Dyson's  equation  for  the  surface  Green's 
function  takes  the  following  form: 


G(m,n)  = G°(m-n)  +G°(m+1)  V (-1 ,0)G  (0  ,n)-*G°  (m+1  )V  (-1 , -1  )G  (-1 , n) 


-*G°  (m)V(0,-l)G  (-1  ,n) -tG°  (m)V  (0 ,0)G  (0  ,n) 


(A.  3) 


In  writing  down  the  above  equation  we  used  the  property  that  due  to  translational 
invariance,  the  bulk  Green's  function  G°  depends  on  the  indices  m and  n only 


./  • ' ' n.  1 ,'^1 


through  the  different  m-n.  Taking  into  account  the  fact  that  the  surface 
Green's  function  G(m,n)  vanishes  if  m and  n refer  to  planes  on  opposite 
sides  of  the  interface,  we  can  solve  Eq.  (A. 3)  for  the  general  matrix  elements 
of  the  surface  Green's  function.  This  is  given  by: 


G (m,n)  = G°(m-n)  + [g°  (m+l)V(->  ,0)  + G°  (m)V(0,0)]  G°(-n)  (l-G°(l)  x 


V(-1,0)-G  (0)V(0 


,o)]"1 


(A.  4a) 


for  m,n>0 


G(m,n)=G°(m-Ti)+  [g°  (m)V(O, -1)  + G°  (m+1  )V (-1 , -1 ) G°(-l  -n)  [l-GC(-l)  x 


V(0,-1)  -(f(-l)V(-l,-l)|  _1 


(A4.b) 


for  m,n  <-l 


G°(n) 


4 ( 


* \ 2T 


Using  this  expression  and  the  expressions  (A.l)  - (A.  2)  for  the  perturbation, 
we  obtain  from  Eqs.  (A. 4)  that  the  matrix  elements  of  the  surface  Green's  func- 


tion are  given  by: 


G(m,n)  “ e 2T  j 


m+n  iu,+(u-'-2U) 
in-(o)-2U) 


(A.  5a) 


A- 2 


form,n<-l,  and 


G(m,n)  =0  (A- 5c) 

otherwise. 


FIGURE  CAPTIONS 


Figure  1 The  local  density  of  states  of  a (100)  interfaced  crystal.  The 

difference  in  the  Fermi  energies  of  the  two  metals  is  2 E. . (a) 

The  first  atomic  layer,  n = 0.  (b)  The  second  atomic  layer,  n = 1. 

(c)  The  third  atcmic  layer,  n = 2.  The  shifted  bulk  density  of  states 
is  shown  as  a dashed  curve. 

Figure  2 The  local  density  of  states  of  a (100)  interfaced  crystal.  The 

difference  in  the  Fermi  energies  of  the  two  metals  is  SE^.  (a)  The 
first  atomic  layer,  n = 0.  (b)  The  second  atomic  layer,  n = 1. 

(c)  The  third  atonic  layer,  n = 2.  The  shifted  bulk  density  of 
states  is  shown  as  a dashed  curve. 

Figure  3 The  interface  bound  states  energies  as  a function  of  the  coupling 
constant.  The  full  and  the  dashed  curves  correspond  to  AE^  = 0, 
and  AEf  = 2T,  respectively. 

Figure  4 The  change  in  the  integrated  density  of  states  due  to  the  creation  of 
the  interface. 

Figure  5 The  interface  energy  of  transition  metals  having  na  and  d elec- 
trons. (a)  n = 0,1,2,  (b)  n = 3,4,5. 
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